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Abstract 

The description of electron current through a splitting is a math- 
t — ■ ematical problem of electron transport in quantum networks [T]. 

' For quantum networks constructed on the interface of narrow-gap 

semiconductors |30[ E] the relevant scattering problem for the multi- 
dimensional Schodinger equation may be substituted by the corre- 
. sponding problem on a one-dimensional linear graph with proper self- 

adjoint boundary conditions at the nodes [THl El ESI ESI \H\ EQ1 HI ESI 
Ell EH El El [I]- However, realistic boundary conditions for splittings 
have not yet been derived. 

Here we consider some compact domain attached to a few semi-infinite 
lines as a model for a quantum network. An asymptotic formula for the 
scattering matrix for this object is derived in terms of the properties of 
' the compact domain. This allows us to propose designs for devices for 

. manipulating quantum current through a splitting El [HI EH HU El] • 
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Introduction: current manipulation in the resonance 
case 



In this paper we discuss the scattering problem on a compact domain with a 
few semi-infinite wires attached. This is motivated by the design of quantum 
electronic devices for triadic logic. In the papers [3], [15] a special design of the 
one-dimensional graph which permits manipulation of the current through 
an elementary ring-like splitting is suggested. This permits, in principle, 
^ . manipulation of quantum current in the resonance case to form a quantum 

switch. Another device for manipulating quantum current through splittings 
is discussed in [231 E] • I n [22] the special design of the splitting formed as a 
circular domain with four one-dimensional wires attached is used to produce 
a triadic relay. 

In order to illustrate the basic principle of operation consider the self-adjoint 
Schrodinger operator 

C = -A + q(x) 

2*1 - n 
an \an ~ u - 

on some compact domain Q. In this paper we will only consider the case 
O C K 3 |23j (for other cases see also [31 [151 E2])- Roughly speaking the 
solution of the Cauchy problem 

idt ~ ^ (i) 
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is given in terms of eigenfunctions <p n 

¥(x,t) = ^a n e !A "V„(i). 

n 

Picking a specific mode ipo with energy Ao we suppose that ipo disappears 
on some subset Iq C 0. Connecting 'thin channels' at various points on the 



v.. 




figure 1. Resonance switch 



boundary of and introducing an excitation of energy Ao along the chan- 
nels we can hope to create a switching effect. Essentialy this is achieved by 
varying q{x) so that Iq D dO, coincides with the connection point of a 'thin 
channel'. 

Implicit in our construction is the assumption that the energy of the elec- 
trons in the device is equal to some resonance eigenvalue of the Schrodinger 
operator on Q. We refer to this as the resonance cas$\. 
Another assumption which we have made above is that Ao is a simple eigen- 
value of C. We will show that the case of multiple eigenvalues is a simple 
generalisation of the case for simple eigenvalues, see [3| [15]. 
In the first section we give a brief description of the connection of the thin 
channels (here they are modelled by one-dimensional semi-lines) to the com- 
pact domain, for more details see [23]. In the second section we derive an 
asymptotic formula for the scattering matrix in terms of the eigenfunctions 
on the compact domain. In the last section we briefly discuss some simple 
models of a quantum switch constructed on the basis of this asymptotic 
formula. 



This has interesting implications when we consider the effect of decreasing the length 
scale — or equivalently scaling up the energy — viz. the effect of non-zero temperature 
becomes negligable for sufficiently small length scales, see [14] . 
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1 Connection of compact domain to thin channels 



As we mentioned above the thin channels are modelled by one-dimensional 
semi-lines. This is justified by an appropriate choice of materials (narrow- 
band semiconductors) and energies, see pEH El [22]. We assume that these 
channels are attached at the points {a±, a-i, ...,azv} C 30 (perturbation of 
the operator £ at inner points {ajv+i, ajv+2, ■ ■ ■ , Q-n+m} may be considered 
using the same techniques as for {ai, a^, ...,ajv} [2 El E2] although we do 
not consider this here). 

We refer to £, defined above, as the unperturbed Schrodinger operator. £ is 
restricted to the symmetric operator £q defined on the class -Do of smooth 
functions with Neumann boundary conditions which vanish near the points 
a\,a2, ■ ■ ■ , ajv- The deficiency subspaces, Af±i, of the restricted symmetric 
operator £q, 

[£* ±i] 4, = o 

for complex values of the spectral parameter A coincide with Greens func- 
tions G\(x,a s ) of £ which are elements of L2(f2) but do not belong to the 
Sobolev class W^iS^)- In the case when Q is a compact one-dimensional 
manifold (a compact graph) these Greens functions are continuous and can 
be written in terms of a convergent spectral series [3|. However, when 
f2 C R 2 , M 3 the deficiency elements will have singularities and we must use 
an iterated Hilbert identity to regularise the values of the Greens function 
at the poles. 

It is well known, for C M 3 , that the Greens function admitts the repre- 
sentation inside O 

i\f\\x-y\ 

G° x (x,y) = — ■ + g(x,y, A) (2) 

4-7r|x — y\ 

where g(x, y, A) is continuous. The potential-theory approach gives the 
asymptotics of the Green function near the boundary point a s € <9S7 

G° x (x,a s )~— — -+L s (x)+B s (x,X), (3) 

2tt\x — a s \ 

where L s is a logarithmic term depending only on d£l and B s is a bounded 
term containing spectral information [8]. 

In order to choose regularised boundary values we use the following lemma 
[23j (here we assume C > — 1 is semi-bounded from below): 

Lemma 1 For any regular point A from the complement of the spectrum 
o~(£) of £ and any a € {a s }^_^ M the following representation is true: 

G\(x, a) = a) + (A + 1)GU * G x (x, a), 

where the second addend is a continuous function of x and the spectral series 
of it in terms of eigenf unctions (pi of the nonperturbed operator £ 

(A + i) -,. G> (,,.)-(A + i)i: (A , w + ( ff A) 

is absolutely and uniformly convergent in Q. 
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The proof of this lemma is based on the classical Mercer theorem along with 
the Hilbert identity [23] . 

It is well known that the domain of Cq can be written as the direct sum 

D* = Do+Mi+M-i (4) 

so for any u £ Dq 

u = u + ^ J 4+G i (x,a s ) + ) j A~G-i(x,a s ). 

s s 

We define u € Dq in terms of the coordinates 
A s = A^ + A s , 
B s = lim 

x— >a s 

the singular and regular amplitudes respectively since it is clear from the 
above lemma that A s is the coefficient of the singular part and B s the 
coefficient of the regular part of u G Dq. The boundary form of Cq may be 
written in terms of A s B s as a hermitian symplectic form 

(/>, v) - (u, C* v) = £ B u s Al - A U S B V S . (5) 
1.1 Self-adjoint extensions 

We recall that to each boundary point a s , s = 1, . . . , N, there is attached a 
semi-infinite ray. On the s-th ray we define the symmetric operator 

d 2 

h,o = + q s (x s ), 

on functions which vanish at x s = (which is identified with a s £ fl). 
Let us consider the symmetric operator Co ii } o © ^2,0 © ••• © In,o- The con- 
nection between the compact domain and the rays is given by (a particular) 
self-adjoint extension of this operator. The boundary form of the adjoint 
Cq © Z* © Z|,o © ••• © ^n,o is easily seen to be 

N N 

£ {B^A* - A*Bf) + J] K(0)W7(0) - u s (0)VM) ■ (6) 

s s 

It is well known that the self-adjoint extensions of Co © Z^o © h,o © ••• © 
In,o correspond to Lagrange planes in the Hermitian symplectic space of 
boundary values equipped with the above boundary form [27]. In general, 
if A, B are (vectors of) boundary values for some symmetric operator then 
any self-adjoint extension can be described by 

l(U-l)A + ±(U + l)B = Q 



u{x) - ^2 A t UGi(x, a t ) 
t 
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for some unitary matrix U [141 [T3] . 

We choose the particular family of self-adjoint extensions which correspond 
to the following boundary conditions at the N points of contact of the rays 



A s \ = (0 p\( B s 
u 8 (0) J V / V 



(7) 



for s = 1, . . . , N and /3 > 0. The resulting self-adjoint extension we denote 
by Cp. The parameter (3 is a measure of the strength of the connection 
between the rays and the compact domain — in the limit (3 — ► the resolvent 
of jC/3 converges uniformly to the resolvent of C on each compact subset of 
the resolvent set of C [23] . 

2 Asymptotics of the scattering matrix 

For the remainder we assume that the potential on the rays q s {x s ) = is 
zero. 

We use the ansatz 

u s = f s (x s , -k)6 a i + f s (x s , k)S 8 i, (8) 

for the scattered wave generated by the incoming wave from the ray attached 
to the point ai. Here f s (x s ,±k) are the Jost solutions [28J, in this case 
(q s (x s ) = 0) just the exponentials 

fs{ x si ifc) = e s , 

and A = k 2 is the spectral parameter. 

^From the boundary conditions ((ZJ) we get N equations 

A s = pf' s (0,-ik)6si+(3f' s (0,ik)S sl 
f3B s = f s (0,-ik)5 sl +f s (0,ik)S sl . (9) 

Inside f2 the eigenfunction u(x, k) may be written as a sum of Greens func- 
tions at the spectral parameter A = k 2 

N 

U 



(x, k) = y~]C s G\(x,a s ). 



Using the Cayley transform between the spectral points i and A one gets a 
relationship between these Greens functions and the deficiency elements (as 
defined above) so that 



lim [G(x,a s ,\) - $tG(x,a s ,i)] = [ \ + X ^ Gi(a s ),Gi(a s ) ) =g s (X). 



x^a s \L — All 

Consequenttly we can show that u has the following asymptotics as x — > a s 
u - C a StGi(x, a s ) + C s g s (X) + C t G x (a s , a t ) + o(l). (10) 
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It follows that for the scattering wave the symplectic variables are related 
by 



A s 
B s 

that is B = QA where 



C s 

g s (X)C s + ^2C t G x (a s ,a t ), 

t^s 



Putting this into 



/ 9 l W 

Gx(a 2 ,a 1 ) 

\ G\(a N ,ai) 



G\(ai,a 2 ) 
5 2 (A) 



G\(ai,a N ) \ 
G\{a 2 ,a N ) 

9 N W J 



(11) 



we can solve for the scattering matrix to get 

I + ik(3 2 Q 



S 



I-ik(3 2 Q' 



(12) 



Let us choose an eigenvalue Ao of the unperturbed operator C on ft. We 
suppose that Ao has a p-dimensional eigenspace, which we denote IZo, with 
orthonormal basis {y?o,i} p - The following important technical statement 
close to Lemma 1 above is true 



Theorem 1 The elements of the Q-matrix have the following asymptotics 
at the spectral point Aq • 



An - A 



+ Qo{a s ,at,X), 



where Qo{a s , at, A) is a continuous function at the point A = Ao- 

We will use this result to prove an asymptotic formula for the scattering 
matrix in the limit of weak connection between the compact domain and 
the rays. 

Consider the mapping V : L 2 {Q) -> which gives the vector of values 
of a function in L 2 (Q) at the nodes of each of the iV rays. To distinguish 
between functions and elements of we use the notation 



and we denote 



Proposition 1 It is possible to choose an orthonormal basis {^o,i} p forTZo 
which forms an orthogonal, but not necessarily normalised, basis for Rq 
under V . 
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Proof: Given some orthonormal basis {(po,i} p f° r T^o we see that 



'■>0,i 



is also an orthonormal basis where U £ U(p). 
The inner product of the image under P 



r,s=l 



shows that finding an orthogonal basis for TZq amounts to finding the unitary 
matrix U which diagonalises A rs = (<po r \<Po,s) • 



□ 



This allows us to write Q in 'diagonal' form 
1 



Q 



A - A u 
+ Qo(A) 
A 



■O,l)<0O,ll 

+ Qo(A) 



A -A 

where m < p is the dimension of Rq. 



(13) 



Theorem 2 If Xq is an eigenvalue of C then for vanishing coupling (3 ~ 
the scattering matrix of Cp has the form 



5(A ) = -I + 2P -2^(^o/3 2 ^o ± Qo^o ± 

s=l 

= -I + 2P + 0(/? 2 ) 
where Pq is the orthogonal projection onto Rq. 
Proof: Using equation (fl~3l) . 



(14) 



5(A) 



i/c/? 2 A) , „ 2 ~ ' 
1+ , " + ik0 2 Q Q 



An - A 



I : : ikp Q 



-1 



Since Do = the ma trix Eq = I 
denominator can be written 



A„-A 



A -A 

is invertable. Consequently the 



_ ik(3 2 D 2 ' 

I - "T r- - ikp Qo 

An — A 



= [[I-ik^QoE^Eo] 1 

= E^ 1 [l-ik^QoEu 1 ]- 1 ■ 

Again the matrix I — ikj3 2 QqEq 1 has an inverse for A ~ An since Qo = Qq. 
This gives the following expression for the scattering matrix 

S{\) = - [EqEq 1 + ikpPQoEZ 1 ] [I - ik^QoEv 1 ] -1 
= - [EqEq 1 + ik(3 2 Q E^] Y,(ikf3 2 Q E^) s . 

s=0 
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Denoting 7Tj = sj (0o,i|0o,«) and diagonalising we can write, 



i / ifc/3 2 7r? ikd 2 ^ 2 N ' 

,. , Aq — A Aq — A N 

aiag I A - A - ikfj^l' ' " ' A - A- ^/3 2 vr 2 ' A ' ' ' ' ' 1 



' m 



Therefore 



Furthermore 



Urn E^=P^=I-P . (15) 

A— >Ao 



f' A -A + i/c/? 2 7T 2 Aq - A + ikft 2 ^ 
00 _diag i Ao-A-^Trf-'Ao-A-^/? 2 ^' 



lim P*^" 1 = Pq" — Po = I — 2P . (16) 

A— >Ao 



which gives us the limit 

iFvom these limits we get 

S(X ) = -[l-2P + ikof3 2 Q Pj-]Y,(i k oP 2 QoP ± ) s 

s=0 

= —I + 2P — 2 ^(^o/^P^QoPo 1 -)* 

s=l 

= -I + 2P + O(/3 2 ). □ 



This formula appears to imply that there may be non-zero transmission 
in the case of zero connection between the rays. Actually the transmission 
coefficients are not continuous with respect to A uniformly in (3 [3J [23] . The 
physically significant parameters of the system are obtained by averaging as 
functions of A over the Fermi distribution so that there is no transmission 
for /3 = 0. 

Corollary 1 If Xq is an eigenvalue of C such that Pq = I then the above 
formula is independent of (3, ie. 

S(X ) = I 

Consequently, when we have pure reflection at an eigenvalue of the unper- 
turbed operator, we have pure reflection regardless of the strength of the 
interaction between the rays and the compact domain. 



3 Simple models 

In [22] the authors discuss the case where f2 is the unit disc in P 2 and there 
are four one-dimensional wires attached at the points ip = 0,7r,±7r/3. The 
dynamics on f2 is given, using polar coordinates (r, 9), by the dimensionless 
Schrodinger equation 

- + [V + er cos(0)] * = A^ (17) 



8 



on the domain with Neumann boundary conditions at the boundary: 

d% 

dn 

The dimensionless magnitude e of the governing field is choosen so that the 
eigenfunction corresponding to the second smallest eigenvalue has only two 
zeroes on the boundary of the unit circle which divide the circumference in 
the ratio 2:1. It is then easy to see that by rotating the potential V one 
may redirect the quantum current from the wire attached to the point ip = 
to any other wire with all of the other wires blocked [22J. 
The analysis in this case is similar to the analysis given above except there is 
now only a logarithmic singularity in the Greens function and the Krein for- 
mula for infinite deficiency indices |18} l2l] and infinite-dimensional Rouchet 
theorem |12| play a central role. A large amount of the calculation was done 
using Mathematica. 

In [TS], using the above asymptotic formula for the scattering matrix 
to choose appropriate parameters, the author discusses the case where VL 
is simply a one-dimensional ring and there is an angle of ir/2 between the 
rays — see figure 2. a). By applying a uniform field to the ring, q = for the 
open state and q = — 3 for the closed state, it is easy to see that a switching 
effect is produced where the Fermi energy is assumed to correspond to the 
smallest eigenvalue of the unperturbed operator on the ring, ie. Ao = 1. See 
also [7] where a similar construction is considered. 
Another possibility is to consider a device — see figure 2.b), the angle be- 




a) Interference b) Potential barrier 

switch switch 

figure 2. 



tween the rays is now tt — with similar parameters where now we switch the 
current by raising a potential barrier, q = for the open state and q = 3 
for the closed state, instead of using interference effects. Clearly, unlike the 
first case, the efficiency of such a switch will be limited by tunneling. 
A more detailed discussion of the properties of these two models (including 
plots of the averaged conductance in the closed, a c , and open, a Q , states) is 



= 0. 

r=1 



9 



presented in [16]. Here we just present the main observations, viz. for both 
switches the open state — possibly due to tunneling effects — is more difficult 
to achieve; and in the limit of small (3, the properties of the switches im- 
prove. This is probably due to the fact that weak coupling between the ring 
and rays improves the open state of the switches. On the other hand, in 
the limit (3 — »■ 1, the ratio d c ja for the second example rapidly decreases 
to a bound due to tunneling which may be calculated from the transmission 
coefficient 

« 4.57 x 10 3 . 

/3=1 

The first switch does not have this bound and consequently for sufficiently 
low temperature or small radius (see first footnote) we conjecture that it 
will have better properties. 



lim — 
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